On the turbulent a-disks and the intermittent activity in AGN 
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ABSTRACT 

We consider effects of the MHD turbulence on the viscosity during the evolution of the 
thermal-viscous ionization instability in the standard a-accretion disks. We consider the pos- 
sibility that the accretion onto a supermassive black hole proceeds through an outer standard 
accretion disk and inner, radiatively inefficient and advection dominated flow. In this scenario 
we follow the time evolution of the accretion disk in which the viscosity parameter a is con- 
stant throughout the whole instability cycle, as implied by the strength of MHD turbulence. 
We conclude that the hydrogen ionization instability is a promising mechanism to explain the 
intermittent activity in AGN. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: evolution 



1. Introduction 

The standard accretion disk (Shakura & Sun- 
yaev 1973) is known to be subject to the thermal - 
viscous instability due to the partial hydrogen ion- 
ization (Meyer & Meyer-Hofmeister 1981; Smak 
1984). As a result of this instability the disk cycles 
between the two states: a hot and mostly ionized 
state with a large local accretion rate and a cold, 
neutral state with a low accretion rate. This insta- 
bility has originally been proposed to explain the 
large amplitude luminosity variations observed in 
cataclysmic variables (CV; Smak 1982; Meyer & 
Meyer-Hofmeister 1982). It is also believed that 
the same mechanism is responsible for the erup- 
tions in soft X-ray transients (SXT; see e.g. Can- 
nizzo, Ghosh & Wheeler 1982; Dubus, Hameury & 
Lasota 2001; Lasota 2001 for review). The ioniza- 



tion instability was also shown to operate in the 
disk around a supermassive black hole in active 
galactic nuclei (AGN; Lin & Shields 1986; Clarke 
1988; Mineshige & Shields 1990; Siemiginowska, 
Czerny & Kostyunin 1996). 

The characteristic timescales of a cycle activity 
scale roughly with the mass of a compact object 
(Hatziminaoglou et al. 2001). Therefore, the ob- 
served timescale of cycle of order of years in bina- 
ries translates into thousands to millions of years 
in galaxies, which harbor a supermassive black 
hole. 

The variability amplitudes and timescales ob- 
served in Galactic binary systems require the vis- 
cosity to be smaller in the quiescent disk than dur- 
ing the outburst. Various ad-hoc viscosity scaling 
laws have been used in the disk evolution models, 
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including the simplest one with a constant viscos- 
ity, however usually a is assumed to be 4-5 times 
larger during the outburst than in the quiescence 
(Cannizzo 1993). The value of ahot determines 
the timescale of an outburst, while a co id governs 
the separation between the subsequent outbursts. 
The difference between these two parameters de- 
termines the outburst amplitude. 

The dependence of the viscosity parameter on 
the disk state is suggested to come out as a prop- 
erty of the magneto-rotational instability (MRI) 
mechanism (Balbus & Hawley 1991; 1998). Non- 
linear development of this instability, which is as- 
sumed to be a primary source of the disk turbu- 
lence, is sensitive to the presence of resistive diffu- 
sion of magnetic field (Hawley, Gammie & Balbus 
1996), measured by the magnetic Reynolds num- 
ber ReM- Gammie & Menou (1998) calculated the 
Reynolds number for the two states of the disk in 
CVs and show that it is low in the quiescence. It 
means that the magneto-hydrodynamical (MHD) 
turbulence dies away and the matter accumulates 
in the outer disk. 

Menou & Quataert (2001) argued that the 
Reynolds number will not be low enough to sup- 
press the MHD turbulence in the quiescent disk 
of AGN. Therefore the efficiency of angular mo- 
mentum transport should be comparable in the 
hot and cold states (a co id w ahot)- hi such a case 
the outburst amplitude is dramatically reduced, 
and the ionization instability leads only to a small 
amplitude flickering (Sicmiginowska et al. 1996). 

The conclusion about the small amplitude of 
outbursts was obtained, however, for a geometri- 
cally thin, optically thick disk extending down to 
the marginally stable orbit. This assumption may 
not be correct. When the local accretion rate is 
low, the local radiative cooling may not be efficient 
(Rees et al. 1982; Bcgclman 1985). The innermost 
part of the disk may be replaced with an optically 
thin, hot and possibly two-temperature plasma. 
An example of such a solution, advection domi- 
nated accretion flow (ADAF) was introduced by 
Ichimaru (1977) and Narayan & Yi (1994). The 
development of the inner flow through the disk 
evaporation was recently discussed by several au- 
thors, starting with Meyer & Meyer-Hofmeistcr 
(1994). 

In this paper we first perform a self-consistency 
check of the a description of the angular momen- 



tum transfer. Having justified a co id = «hot ap- 
proach to AGN accretion disks, we consider the 
time evolution of the disk under the influence of 
the ionization instability. We also take into ac- 
count the evaporation of the inner disk. 

We define a parameter space where the instabil- 
ity zone (i.e. the partial hydrogen ionization zone) 
and the evaporation region overlap. The overlap- 
ping of these two regions strongly enhances the 
outburst amplitude due to the time evolution of 
the disk evaporation radius. We illustrate this be- 
havior with exemplary lightcurves. 

The structure of this article is as follows. In 
Section 2 we define the criterion for efficient de- 
velopment of the MHD turbulence in the accretion 
disk, taking into account its vertical structure. We 
confirm that indeed the bi-modal viscosity behav- 
ior is characteristic for Galactic binary systems. 
However, a single value of the viscosity parameter 
is appropriate in the case of AGN disks. In Section 
3 we check if the zone of partial hydrogen ioniza- 
tion overlaps with the disk evaporation region and 
we calculate the evolution of the disk luminosity 
in case of the evaporated inner disk. Finally, in 
Section 4 we discuss our results and give conclu- 
sions. 

2. Magneto-rotational turbulence 

In this Section we estimate the strength of the 
MHD turbulence in the a accretion disks and in- 
vestigate whether the viscosity prescription with 
a constant parameter a is appropriate for both 
Galactic systems (CVs and SXTs) and AGN. We 
only consider the instability due to partially ion- 
ized hydrogen and do not include regions affected 
by the radiation pressure instability (see Janiuk, 
Czerny & Siemiginowska 2002 for that model de- 
tails). We assume that the viscous stress tensor is 
proportional to the gas pressure P gas and therefore 
there is no radiation pressure dominated branch on 
the stability figures (so-called S-curves) described 
below. 

We use the dimensionless accretion rate in Ed- 
dington units, m = Af/MEdd, assuming the effi- 
ciency of 1/12, as implied for the Schwarzschild 
black hole by the Newtonian potential: 

M Edd = 3.52-^— [M e /yr\. (1) 



2 



The disk vertical structure (indicated by z) is 
calculated by solving the equations of viscous en- 
ergy dissipation, hydrostatic equilibrium, and en- 
ergy transfer: 



dF 
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= aP gas (- 
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where we solve for temperature, density and pres- 
sure profiles, T(z), p(z) and P(z). Here is the 
Keplerian angular velocity, a and c are physical 
constants. F\ is the energy flux transported lo- 
cally in the direction perpendicular to the equa- 
torial plane and carried either by radiation or by 
convection: 

F\ = F ra d V ra d < V a d, F\ — Frad+^conv V ra d > V a d 

(5) 

The frequency-averaged opacity k is taken to be 
the Rosseland mean and includes the electron scat- 
tering, free- free and bound- free transitions. The 
opacity tables are from Alexander, Johnson & 
Rypma (1983) and Seaton et al. (1994). The 
presence of dust and molecules is included in the 
opacity description. The details of the model were 
discussed in Pojmahski (1986) and Rozanska et al. 
(1999). 

Figures 1 and 2 show, in X vs. T c ff (surface 
density vs. effective temperature) plane, the local 
solutions of the disk vertical structure calculated 
for a range of accretion rates for two extreme cases 
of the mass of the central object: 1M Q and 3 x 
10 9 M Q . The solutions located on both lower (be- 
low point A) and upper (above point B) branches 
of the S-curves are stable. The slope of the middle 
branch is, at least in some parts, negative which 
means that at that range of accretion rates the 
disk is thermally and viscously unstable. 

Note, that in case of CVs (Figure 1) the lo- 
cation of the starting point of the instability is 
very sensitive to the adopted value of a. For a 
small viscosity parameter, a — 0.02, the point A 
is at much lower temperature than in case of large 
a = 0.1. The latter 'smoothens' the S-curvc in its 
lower part, so that only a single 'wiggle' remains. 
Therefore the point A, defined as the first critical 



point above which the S-curve slope becomes neg- 
ative, is shifted upwards in temperature in com- 
parison with that characteristic for small viscosi- 
ties. This fact should be taken into account when 
modeling the limit cycle in dwarf novae by means 
of the 'combined' S-curves with a co u < oihot- 

This behavior is however characteristic only for 
accretion disks around 1M Q , while vanishes al- 
ready for SXTs. Also for supermassive black holes 
there are always two 'wiggles' in the S-curve, re- 
gardless of the viscosity parameter. 

The ionization instability is characteristic for 
all accretion disks but the range of accretion rates 
(i.e. the location of the S-curve on the T— £ plane) 
depends on the chosen disk radius. We can invert 
this problem and say that for a fixed external ac- 
cretion rates the instability will appear only for a 
certain range of radii. At larger radii the disk is on 
the lower stable branch while at lower radii it is on 
the higher, also stable branch. Simple analytical 
formulae for the unstable zone were provided by 
Siemiginowska et al. (1996) and applied by Menou 
& Quataert (2001) in their analysis. 

In Figure 3 we show the radial extension of the 
unstable zone as a function of the accretion rate 
(in Eddington units) calculated numerically from 
our disk model. The position of the instability 
zone depends on the mass of the central object so 
we choose values representative for all objects from 
CVs to AGN with a large black hole mass. The 
plot shows the radial extension as a function of the 
accretion rate (in Eddington units) . The horizon- 
tal slices of the shaded region correspond to the 
accretion rates, for which the disk is unstable at a 
particular radius. For m « 0.01 (an accretion rate 
typical for many objects) the instability zone is lo- 
cated at ~ 10 5 i?schw for a CV, at — 2 x 10 3 i?schw 
for M ~ 10 6 M Q and at ~ 500i? S chw for an ex- 
tremely massive black hole of M — 3 x 10 9 M Q . 

Time evolution of the unstable part of the disk 
proceeds roughly in the form of oscillations be- 
tween the upper and the lower stable branches. 
As argued by Gammie & Menou (1998), high tem- 
perature upper branch conditions are always fa- 
vorable for the development of the efficient MRI 
and high values of viscosity are appropriate there. 
Lower branch conditions are different and the MRI 
mechanism may not be efficient. Therefore, we 
perfom our self-consistency check paying particu- 
lar attention to the lower branch solutions at the 
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vicinity of the turning point A. 

2.1. Coupling of the magnetic field with 
the gas 

The behavior of the magnetic field is governed 
by the fluid conductivity a. The time-dependent 
magnetic field in the disk is described by the equa- 
tion (Parker 1979): 

pi ~D 

— = Vx(vxB) + r]V 2 B (6) 
at 

where r] = c 2 /4na denotes resistivity. The char- 
acteristic diffusion time in which the initial con- 
figuration of the magnetic field will decay is equal 
to t — L 2 /t], where L indicates the characteristic 
spatial scale. For timescales short in comparison 
with r the second term in Equation 6 can be ne- 
glected and the magnetic field lines are frozen into 
the gas. Magnetic Reynolds number, defined as: 



Re 



M 



VT 

T 



(7) 



can be used to distinguish between two cases: (1) 
the diffusion of field lines within the disk and 
(2) the lines are frozen in and carried along with 
the matter. Here we identify the velocity v with 
the sound speed in the disk and the length L of 
the magnetic field spatial variations with the disk 
thickness, as commonly used in the simulations 
(see e.g. Hawley et al. 1996; Gammie & Menou 
1998). 

Apart from the ohmic diffusion also the am- 
bipolar diffusion may be important. The ambipo- 
lar Reynolds number is defined as: 



Re A = 



n 



(8) 



where v n \ is the frequency of collisions between 
ions and neutral particles. 

We can calculate the both Reynolds numbers 
locally in the disk. In order to determine the re- 
sistivity, r\ = (c 2 TO c /47re 2 ) x (u en /n e ), we have 
to estimate the number density of electrons n c 
and neutral particles n n - The frequency of colli- 
sions between electron and neutral particles is de- 
fined as v cn = 8.3 x lO-^T 1 / 2 ^ hi s" 1 (Drainc, 
Roberge & Dalgarno 1983). The ionization frac- 
tion x c — n c /n n is calculated from the Saha equa- 
tion for any given value of temperature and den- 
sity in the disk. In this temperature range the 



free electrons originate mostly from the ionization 
of hydrogen, sodium, potassium and calcium. 

The magneto-rotational instability acts over 
the entire disk thickness. Therefore, we need to 
consider the entire vertical disk structure when 
we estimate the coupling of the magnetic field to 
the gas. This is why instead of calculating ReA 
and ReM at the equatorial plane, we introduce 
the density weighted values of these parameters, 
integrated over the disk height: 



j" Re A p(z)dz 



Ret 



So P dz 



and 



Re M p{z)dz 
f" pdz 



Re 



M 



(9) 



(10) 



where p(z) is the density profile determined from 
the vertical structure model and H = z max is the 
disk thickness. 

Numerical simulations (Hawley ct al. 1996; 
Hawley & Stone 1998) show that below the critical 
Reynolds numbers (defined as in Equations 7 and 
8): 

Reft* = 10 4 (11) 

and 

Ref* = 100 (12) 

the action of MHD turbulence is inefficient and 
cannot be important in the angular momentum 
transport in the disk. As a consequence, the vis- 
cosity parameter (a « a m ag) is very small. 

If the disk Reynolds numbers drop below crit- 
ical values in the quiescence (e.g. on the lower 
branch of the S-curve), while rise above the crit- 
ical values in the outburst (e.g. on the upper 
branch of the S-curve) then the assumption that 
c*coid <C ahot (see Figure 1) is justified. Other- 
wise, the viscosity parameter should be the same 
on both cold and hot branches of the S-curve (see 
Figure 2). Below we calculate the Reynolds num- 
bers in the quiescent disk and test them against 
the critical values in different types of accreting 
systems. 

2.2. Results 

2.2.1. Accretion disks in binary systems 

We consider an accretion disk in CVs, e.g. 
around a white dwarf with a mass M=1M Q and 
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calculate the disk vertical structure at two exem- 
plary radii: 10 4 Rs c hw and 10 5 Rs c hw Figure 4 
(a) shows the magnetic and ambipolar Reynolds 
numbers for the case of small viscosity parame- 
ter a = 0.02. It is plotted as a function of the 
accretion rate corresponding to the lower branch 
of the S-curve (e.g. in Figure 1 below the point 
A). Point A indicates the highest accretion rate 
on the lower, stable branch of the S-curve. Both 
Reynolds numbers are below the critical values for 
most part of the lower branch. Magnetic Reynolds 
number is somewhat above the critical value close 
to the point A. The ambipolar Reynolds number 
is well below the critical value at the smallest ra- 
dius and approaches it at the point A only for the 
largest radius. In general, the magnetic number 
at the lower branch is almost independent from 
the radius. The ambipolar number is increasing 
with the radius. However, the value of 10 5 Rs c hw 
is rather an upper limit for the disk size in CV sys- 
tems so any values of the Reynolds number larger 
than those in Figure 4 are not expected. 

In Figure 4 (b) we plot the magnetic and am- 
bipolar Reynolds numbers up to the point A for 
large value of a = 0.1. In this case the Reynolds 
numbers are mostly small, but as the point A 
reaches this time higher temperatures, also their 
values can exceed critical ones in the end of the 
extended stable branch. 

Therefore, whenever the viscosity drops at the 
lower branch it will have a tendency to remain low. 
This result supports the view that the angular mo- 
mentum transport due to MHD turbulence cannot 
be efficient when the disk is in the quiescence and 
the assumption that a co id "C ahot is justified. 

We obtain a similar result for black hole X-ray 
transient systems with Mbh ~ 10M Q as shown in 
Figure 5. The calculated magnetic and ambipolar 
Reynolds numbers fall generally below the criti- 
cal values for the quiescent disk in these systems, 
regarding also the fact that this time there is no 
difference in the location of the point A caused by 
the change in a. 

2.2.2. AGN accretion disks 

Now we consider accretion disks in AGN, e.g. 
systems harboring a supermassive black hole. 
Note, that for an accretion disk around a super- 
massive black hole the temperatures are lower 



than in the stellar mass black hole case. There- 
fore the ionization instability zone is located much 
closer to the central black hole and in our calcu- 
lations we have to consider radii of order of a few 
hundreds i?schw instead of hundreds of thousands. 

We calculate the Reynolds numbers for a range 
of masses between 10 6 M Q and 3 x 10 9 M Q . We 
found that the Reynolds numbers are orders of 
magnitude above the critical values, for the most 
part of the lower branch, and the greater the mass 
of a black hole, the larger Re. In Figure 6 we plot 
magnetic and ambipolar Reynolds numbers as a 
function of the accretion rate corresponding to the 
lower, cold branch of the S-curve for 3 x 1O 9 M . 

Both Reynolds numbers clearly exceed the crit- 
ical values: 10 2 and 10 4 for Re\ and ReM respec- 
tively. Therefore only for very low accretion rates, 
m < 10 -4 , MHD turbulence in the disk would not 
develop and practically the entire cold branch of 
the S-curve is MHD turbulent. This result only 
weakly depends on the adopted viscosity param- 
eter value. For a — 0.02 both Reynolds numbers 
had similar values to these presented in Figure 6 
and were always greater than critical (ReM > 10 6 
and Re a > 10 5 at the disk radius of logi? =16.5 
[cm]) . 

In order to understand whether the above re- 
sult is valid for the entire disk we calculate the 
Reynolds numbers at different locations in the fol- 
lowing way. At each radius we compute the en- 
tire S-curve, such as that in Figure 2, determine 
the position of the turning point A and calcu- 
late the Reynolds numbers ReA(A) and Reu(A) 
at this point. Their values measured along the 
lower branch of the S-curve are lower (equal) to 
ReA{A) and Reyi{A). We summarize this result 
in Figure 7, plotting the Reynolds numbers against 
the radius. Both ReA and ReM increase with ra- 
dius implying that MRI is stronger further out in 
the disk. They decrease towards smaller radii but 
do not drop below the critical values even at 10 

-Rschw 

Our results support the view presented in 
Menou & Quataert (2001) and imply that in the 
thermally unstable AGN disks the viscosity should 
not change between the cold and hot states of the 
disk (i.e. a co id = ahot)- 
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3. Evaporation of the disk 

We concluded above that the viscosity parame- 
ter in AGN accretion disks does not vary between 
the cold and hot states of the disk. As a result the 
ionization instability will cause only small ampli- 
tude luminosity fluctuations not really detectable 
in the observations (Sicmiginowska, Czerny & 
Kostyunin 1996; Menou & Quataert 2001). How- 
ever, variations in the local accretion rate induced 
by the instability may result in the disk evapora- 
tion and affect the disk luminosity (Sicmiginowska 
1998). Here, we investigate the effects of the disk 
evaporation on the overall evolution of thermally 
unstable disk. We assume that a co id = «hot = 0.1. 

3.1. Evaporation radius 

The accretion flow in low luminosity systems 
can occur either via a standard cold disk or via 
an advection dominated flow (ADAF). In the lat- 
ter the bulk of gravitational energy is advected 
into the central object instead of being radiated 
(Ichimaru 1977; Narayan & Yi 1994). The ADAF 
solution forms below a certain transition radius 
that depends on the accretion rate. In classical 
ADAF model (Abramowicz et al. 1995; Honma 
1996; Kato & Nakamura 1998) this relation is: 

-Rcvap = 1.9m~ 2 ao.i-Rschw (13) 

The transition radius between the outer disk and 
inner ADAF was discussed in several papers (e.g. 
Liu et al. 1999; Rozahska & Czerny 2000; Menou 
ct al. 2000; Liu et al. 2002; Meyer & Meyer- 
Hofmeister 2002), under various assumptions. 
Here we adopt the canonical approach. 

The ionization instability can develop only 
when the standard thin disk reaches the tem- 
peratures sufficient for the hydrogen ionization. 
If accretion rates are low the disk can evaporate 
without being ionized (e.g. without reaching point 
A on the stability S-curve) and the instability does 
not develop. Therefore the evaporation radius 
cannot be very large for the ionization to occur: 

-^evap ^ -^instab- (1^) 

In Figure 3 we showed the radial extension of 
the unstable ionization zone as a function of the 
accretion rate. The instability cycle operates be- 
tween the points A and B of the S-curve (cf. Fig- 
ure 2), so we consider this entire range of accretion 



rates to be unstable. In Figure 3 we also indicate a 
location of i? e va P for the ADAF prescription (note, 
that the evaporation radius does not depend on 
mass of a black hole). If the disk evaporates, only 
the unstable region above the solid line may exist 
and contribute to the total disk luminosity. The 
extension of the unstable region does not depend 
on a very much, while the formula defining the 
evaporation radius strongly depends on it. There- 
fore, we should emphasize that our result may be 
very sensitive to the adopted value of the viscosity 
parameter. Possibly, depending on a, the entire 
unstable region of the disk evaporates, or, on the 
other hand, the evaporation radius is very small 
and does not affect the instability. 

Because the location of the evaporation radius 
depends on viscosity, if it can be determined ob- 
servationally it may provide a new method of esti- 
mating the value of the a parameter in accreting 
objects. For example, recently Nayakshin & Sun- 
yaev (2003) proposed a new model explaining the 
X-ray flares in the Galactic Center by the star-disk 
interactions. This model would favor the case of 
the disk extending quite close to the central black 
hole. However, the presence of the disk in the 
Galactic Center is still unclear (see Quataert 2003 
for a review). 

3.2. Luminosity of the disk with inner 
ADAF 

Below we study the evolution of the disk af- 
fected by the ionization instability assuming that 
the viscosity parameter is constant, as favored 
for accreting supermassive black holes. We used 
the method described in Sicmiginowska, Czerny & 
Kostynuin (1996) and consider the two cases: (i) 
the accretion disk extends down to the marginally 
stable orbit and (ii) the central part of the disk is 
evaporated and forms an ADAF, while the stan- 
dard thick disk remains at the outer radii. Since 
the accretion efficiency of an ADAF is low, we as- 
sume that this part of the flow does not contribute 
to the luminosity. 

In order to determine the disk evaporation ra- 
dius during the evolution we compare the local 
accretion rate with the critical (evaporation) ac- 
cretion rate. We use the classical ADAF prescrip- 
tion, expressed by Equation 13. In Figure 8 we 
show radial profiles of the local accretion rate in 
several snapshots during one cycle of the disk in- 
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stability. The straight solid line marks the critical 
accretion rate. Whenever the local accretion rate 
in the disk drops below this critical value, the disk 
evaporates and ADAF is formed. Hereafter, in 
our time dependent calculations we assume that 
the ADAF transition radius is the outermost ra- 
dius at which the disk evaporates. Therefore even 
though the local accretion rate may be equal to 
the critical value at several locations in the disk, 
only the outermost one determines the transition 
radius and the standard disk cannot rebuilt itself 
below this radius. 

In Figure 9 we present lightcurves for the two 
models: (i) the standard disk and (ii) the disk 
with inner ADAF. When the disk extends all the 
way to the marginally stable orbit, as considered 
in the previous time-dependent calculations, the 
luminosity variations are small and only a small 
flickering (Alogidisk ~ 0.2) can be seen in the re- 
sulting lightcurve. This is indicated by the upper 
solid line in Figure 9. 

However, the situation changes dramatically, 
when the evaporation of the inner disk is taken 
into account. In this case the ionization instability 
leads to strong, large-amplitude outbursts. These 
outbursts are due to the large variations of the in- 
ner radius of the accretion disk (e.g. the size of 
the ADAF), as determined by the critical value of 
the local accretion rate. The outburst peak corre- 
sponds to the moment, when the inner disk radius 
is small and the entire accretion disk contributes 
to the total luminosity. The quiescent low lumi- 
nosity state, on the other hand, is achieved when 
the local accretion rate drops below the critical 
value at the large distance from the center and 
most of the inner disk is evaporated. 

4. Discussion and conclusions 

We studied the effects of MHD turbulence and 
the efficiency of angular momentum transport in 
the quiescent accretion disks. We analyzed the full 
disk vertical structure, which enabled us to deter- 
mine accurately the coupling of the accreting gas 
to the magnetic field over the entire disk height. 

Confirming the previous results of Menou & 
Quatacrt (2001), we find that only in the case 
of accretion disks in binary systems, e.g. around 
white dwarfs (CVs) or stellar mass black holes 
(SXTs), the timescale of the magnetic field de- 



cay is too short for the field to be frozen into the 
gas. Therefore, MHD turbulence does not oper- 
ate efficiently and the large amplitude outbursts 
observed in these sources can be accounted for by 
the thermal ionization instability cycle in which 
a co id <C CKhot- On the other hand, in the case of 
accreting supermassive black holes in AGN the gas 
in the disk is well coupled to the magnetic field. 

This effect (ReM ^> 1 and Re a ~3> 1) can be 
only slightly reduced by encounters of free elec- 
trons with charged dust particles. In principle, the 
presence of dust in the disk can affect the number 
density of free electrons. However, the grains evap- 
orate if the temperature exceeds ~ 1800 — 2000 K, 
which is the case near the disk equatorial plane. 
As a result, the role of the dust in suppressing 
MHD turbulence is limited only to the uppermost 
disk layers and cannot have a global effect. We 
conclude, that the AGN accretion disks in quies- 
cence have comparable angular momentum trans- 
port efficiency to the disks in outburst. 

Therefore, we studied the ionization instabil- 
ity assuming a constant value of the viscosity pa- 
rameter a co id = c*hot- We took into account the 
evaporation of the innermost parts of the disk dur- 
ing the quiescence in AGN accretion disks. The 
transition to an ADAF below the evaporation ra- 
dius during the evolution of the instability results 
in large luminosity variations (AlogL^isk ~ 2). 
Because only the outer parts of the disk (R > 
Revap ~ 10 3 i?schw) contribute to the total disk 
emission the minimum luminosity obtained within 
the instability cycle is lower than in the standard 
disk with no evaporation. We note, that also dur- 
ing the outburst the evaporation is important and 
ADAF remains in the innermost regions of the disk 
(Revap ~ 25i?schw), even at quite a high global 
accretion rate, M = O.lMEdd- The evaporation 
radius, however, strongly depends on the adopted 
value of the viscosity parameter a (cf. Equation 
13). 

The a-disks in AGN are self-gravitating beyond 
a certain radius (Paczyhski 1978; Hurc 1998). 
A typical distance of the self-gravitating zone is 
about a thousand Schwarzschild radii (Hure 2000). 
However, the greater the black hole mass and 
lower the parameter a, the closer to the center this 
zone is located. We can include a self-gravitation 
term equal to — AttGY, in Equation 3 of the disk 
vertical structure. This additional term modi- 
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fies the S-curve, so there is a single 'wiggle' on 
the middle branch of the curve, instead of two. 
The middle branch is also broader and reaches 
lower surface densities. The instability happens 
for a narrower range of effective temperatures, 
while the more strongly curved unstable middle 
branch covers broader range of surface density val- 
ues. Therefore, the amplitude of the outburst is 
not expected to be significantly modified in our 
model by the presence of self-gravity. However, 
whether the model actually can be extended into 
the self-gravity region, remains an open question. 
In particular, our description of convection, which 
influences the shape of the S-curve, may not be 
applicable in case of a self-gravitating disk. Also, 
MRI was not extensively studied in the self grav- 
ity region (but see Fromang et al. 2003) and the 
use of viscosity parameter a may not be appro- 
priate here. Instead, self-gravity may either pro- 
vide a mechanism of angular momentum transfer 
(e.g. Paczyfiski 1978), or lead to the disk fragmen- 
tation and star formation (Collin & Zahn 1999). 
The physics of those processes is poorly known. 
Therefore the use of simple models like a-disks 
and comparison of their predictions with observed 
behavior of AGN seems to be justified. 

The variability amplitudes predicted by our 
model are about two orders of magnitude and 
support the observational evidence of intermittent 
activity of AGN. This kind of activity has been 
proposed to explain the formation of giant radio 
galaxies (Subrahmanyan et al. 1996) and the ex- 
tended radio structures in the Giga-Hertz Peaked 
Spectrum (GPS) radio sources (Baum et al. 1990; 
Schoenmakers et al. 1998). The timescales of 
episodic activity in quasars were recently esti- 
mated by Martini & Schneider (2003). Recent 
observations of Compact Steep Spectrum (CSS) 
sources indicate that the evolutionary path of ra- 
dio loud AGN is consistent with an on-off activity 
with duration governed by the mass of the cen- 
tral black hole (Marecki et al. 2003). Also the 
black hole masses implied by the HST observa- 
tions of normal galaxies and AGN suggest that 
every galaxy must undergo an active phase in its 
life (Ferrarese et al. 2001). 

The estimated timescales of the activity cycle 
are within 10 5 to 10 7 years (e.g. Reynolds & 
Begelman 1997 for GPS sources). Di Matteo et 
al. (2003) proposed a feedback mechanism to ex- 



plain the cyclic activity in M87, which occurs on a 
timescale of 10 7 years. In our time-dependent cal- 
culations, we adopted a simplified approach to the 
disk evaporation process, and therefore underesti- 
mated the repetition timescale. This is because 
we do not consider the time of the disk recon- 
struction after it has been evaporated, which in- 
creases the viscous timescale by a factor of ~ 2. 
Assuming that the outburst starts at the outer 
edge of the instability zone we obtain the viscous 
timescale of i visc = l/(aO) x (r/H) 2 w 0.9 Myrs 
for M B h = 1O 8 M . Thus the outbursts will be 
less frequent than these presented in Figure 9. 
However, because during the outburst there is no 
ADAF and the entire disk is present, the decay 
timescale is calculated properly and the duration 
of an outburst is correct. On the other hand, the 
outburst profiles will be less steep than the ones 
calculated here, since the disk density rises grad- 
ually at the beginning of the instability cycle. Al- 
though these two effects are important, the ex- 
tension of the unstable region is not very large in 
comparison with the disk size and the uncertainty 
involved in our calculations is not crucial. The 
detailed model of the time evolution of the accre- 
tion disk with transition to ADAF is clearly worth 
investigating and we address it to the future work. 
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Fig. 1. — The stability curves calculated for the 
disk around M = 1M© at the radius R = 
10 5 i?schw The viscosity parameters are a = 0.1 
(left curve) and a — 0.02 (right curve). The dashed 
line marks the combined stability curve, which 
would result from the assumption that a co id = 
0.02 (on the lower stable branch) and «hot = 0.1 
(on the upper stable branch) . Point A indicates the 
lowest location of the unstable disk and a starting 
point of the ionization instability cycle and the 
instability ends in the point B. 
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Fig. 2. — The stability curve calculated for the disk 
around a supermassive black hole of M = 3 X 
10 9 M Q at the radius R = 300i?schw The viscosity 
parameter is a = 0.1. The point A is the starting 
point of the ionization instability cycle and the 
instability ends in the point B. 
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Fig. 3. — The radial extension of the ionization in- 
stability zone as a function of the accretion rate. 
The contours correspond to the turning points A 
and B on the S-curve (see Figure 2). The thick 
solid line marks the transition radius, resulting 
from the ADAF prescription (see text). The black 
hole mass is M = IMq (top contour), M = 
10 6 M Q (middle contour) and M = 3 X 1O 9 M 
(bottom contour). The viscosity parameter is a = 
0.1. 
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Fig. 4. — The density weighted magnetic and am- 
bipolar Reynolds numbers calculated for the cen- 
tral mass M = IMq, plotted for two disk radii 
R = 10 4 i?schw (dashed line) and R = 10 5 i?s c hw 
(solid line), as a function of the accretion rate 
in the range corresponding to the cold branch of 
the S-curve (cf. Fig. 1). The viscosity is a = 0.02 
(panel a) and a = 0.1 (panel b). 
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Fig. 5. — The density weighted magnetic and am- 
bipolar Reynolds numbers calculated for black 
hole mass M = 10M Q , plotted for two disk radii 
R = 10 4 i?s c hw (dashed line) and R = 10 5 i?schw 
(solid line), as a function of the accretion rate in 
the range corresponding to the cold branch of the 
S-curve. The viscosity is a = 0.1. 
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Fig. 6. — The density weighted magnetic and am- 
bipolar Reynolds numbers calculated for black 
hole mass M = 3 x 10 9 Af Q , plotted for two 
disk radii R — 300i?schw (dashed line) and R = 
500i?schw (solid line), as a function of the accre- 
tion rate in the range corresponding to the cold 
branch of the S-curve (cf. Fig. 2). The viscosity is 
a = 0.1. 
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Fig. 7. — The radial profiles of the magnetic (cir- 
cles) and ambipolar (triangles) Reynolds numbers 
in the turning point A on the S-curve. The points 
mark the values obtained for 10, 50, 100, 200, 300, 
500, 700 and 1000 i?schw at each curve respec- 
tively. The black hole mass is M = 3x 1O 9 M 
and viscosity is a = 0.1. 
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Fig. 8. — The local accretion rate during the cy- 
cle of the disk evolution caused by the ionization 
instability. The thick solid line marks the critical 
accretion rate, below which the disk is evaporated, 
resulting from the ADAF prescription (see Eq. 13). 
The black hole mass is M — 10 8 M© and viscosity 
a = 0.1. 
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Fig. 9. — The disk lightcurve during the evolution 
caused by the ionization instability with constant 
viscosity assumed (a = 0.1). The upper curve 
shows the luminosity flickering resulting from the 
evolution of the whole disk. The lower curve shows 
the sharp outbursts of the disk, which inner part 
is evaporated. The black hole mass is M — 10 8 M Q 
and the global accretion rate is 0.1 MEdd- 
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